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Introduction 



I . 

This report describes position estimation procedures that 
are based on models which relate positional uncertainty to 
various measurement and estimation errors. The procedures were 
developed to be used in analyses that relate the effect of posi- 
tional uncertainty on tactical performance through such factors 
as weapon accuracy as well as to be used operationally. 

In the models, positions are on a plane surface (flat earth 
model). Because of this condition, the models are not intended 
for use in situations where the earth's figure is significant. 

In addition, positional errors are determined by independent 
normally distributed random variables with known means, variances 
and covariances. The support for this condition, other than its 
mathematical convenience, is that it has been used by others, 
for example, see Reference 1. To use the models, one is required 
to specify the means, variances and covariances of the error 
random variables. For these models, this can be done by speci- 
fying the systematic errors (biases) and the error magnitudes 
(standard deviations) . 

The procedure that is described in Section II relates position 
estimates that are based on bearings on or from stations to 
station position uncertainty. One application for the model is 
the analysis of the effect of sonobuoy position uncertainty on 
position estimates determined with passive directional sonobuoys . 

The procedure that is described in Section III relates posi- 
tion estimates that are based on lines of position to line of 
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position uncertainty. In addition to estimating the uncertainty 
in position estimates based on celestial observations, the pro- 
cedure could be used to determine error ellipses for LORAN fixes 
if the standard deviation values required by the procedure could 
be obtained. 

The procedure that is described in Section IV can be used 
to combine position estimates from various sources. The proce- 
dure which is based on conditions that should not be too restric- 
tive in most cases provides both a composite position estimate 
and error ellipse. 
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Station Position Uncertainty and Position Estimates 



II . 



The procedure that is described in this section relates 
position estimates that are based on bearings on or from stations 
to station position uncertainty. The procedure is based on a 
model that is an extension of one that is described on Reference 
2. The model is defined as follows: Each station position error 

is determined by an independent bivariate normal distribution 
with a zero mean vector and a known covariance matrix. Observed 
bearing lines on or from a station are parallel to true bearing 
lines. The distance of each observed bearing line from its true 
bearing line is determined by an independent normal distribution 
with a zero mean and a standard deviation o. 

Because a station's position error is determined by a bivar- 
iate normal distribution, the perpendicular distance s between a 
line at the assumed location of an observed bearing line and the 
observed bearing line is determined by a normal distribution. 

The relation between the bivariate normal distribution that 
describes the station's position uncertainty is indicated in 
Figure 1. In the figure, the positive y-axis direction is north, 
the positive x-axis direction is east, and the origin of the 
coordinate system is at the assumed station position. The x'y'- 
coordinate system is oriented so that the positive x'axis is 
coincident with the major axes of the elliptical contours of 
the bivariate normal distribution that determines the station 
position error and so that the bearing 6 of the positive x 1 -axis 
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Figure 1. The geometry associated with the determination of a s , 
the standard deviation of the normal distribution 
that determines the distance between the assumed 
location of an observed bearing line and the observed 
bearing line. The assumed (mean) position of the 
station is at the origin. The ellipse represents a 
contour on the probability density surface of the 
bivariate normal distribution that determines the 
station's position. 
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satisfies the condition: 0° < 5 180°. The anqle is an 

observed bearing. The x"y "-coordinate system is oriented so 
that the positive y"-axis is in the direction of the observed 
bearing line. As a consequence of these relationships. 



a 2 = a 2 sin 2 ( 6 - 0 ) + 0 2 cos 2 (6- e) 

S X y 

2 2 

where o^, and 0 ^, are the elements of the station position 
error covariance matrix relative to the x ' y ' -coordinate system. 

A procedure for determining position estimates that are 
based on bearings on or from stations at known positions is 
described in Appendix 1 of Reference 3. The procedure is based 
on a model that is equivalent to one that is described in Ref- 
erence 2. The procedure that is described in this section is 
based on a model that is an extension of it. 

The procedure in Appendix 1 of Reference 3 is based on a 
model in which a station's bearing error is determined by a 
normal distribution with zero mean (bias) and standard devia- 
tion e. The bearing error is related to the distance on a cir- 
cular arc between a station's true bearing line and observed 
bearing line. The arc is on the circle with its center at the 
station that passes through an initial estimate of an object's 
position. This distance is determined by a normal distribution 
with mean zero and standard deviation a = re where r is the 
range of the initial estimate from the station and the standard 
deviation e is measured in radians. In the model, arc distance 
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is approximated using a first order approximation which in 
effect replaces the circle with its tangent line at the initial 
estimate. As a consequence, the distance u on the tangent line 
between the observed bearing line and the true bearing line is 
determined by a normal random variable with mean zero and stan- 
dard deviation a. The distance u can be expressed in terms of 
w, the distance on the tangent line between the observed bearing 
line and the initial estimate that is also determined by a normal 
random variable with standard deviation a, and v, the distance 
on the tangent line between the true bearing line and the initial 
estimate. And, as shown in Appendix 1 of Reference 3, this 
distance can be expressed in terms of the unknown coordinates 
of the object's position. 

The effect of station position uncertainty is accounted for 
by the distance s between the observed bearing line and the 
assumed observed bearing line. In the model s is determined by 
a normal distribution with mean zero and standard deviation a 

s 

as given above. This approximation is consistent with the first 
order approximation of arc distance. As a consequence of these 
two approximations, all of the bearing lines are replaced by 
lines parallel to the line joining the initial estimate's 
position and the station's assumed position, both of which are 
known positions. The geometry involved in shown in Figure 2. 

The modified relationships resulting from the introduction of 
station position uncertainty are shown in Figure 3. The 
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Figure 2. The geometry of the position estimation model. The 
bearing 3 and range r of the initial estimate from 
the assumed station position have the role of S 
and r in Appendix 1 of Reference 3. 
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assumed location 
of the observed bearing line 




Figure 3. The quantities u, v and w correspond to the quantities 
u, v and w in Appendix 1 of Reference 3. The auxiliary 
quantities u' and w' are defined in the figure. The 
replacement of by 0“ + a s ^ in the procedure in 
Appendix 1 of Reference 3 is justified by noting that 
with this substitution u' and w' are equivalent to u 
and w. 
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important thing to note is that u' ^ u - s is determined by a 

normal distribution with mean zero and standard deviation 
2 2 t 

[o + o s ] 2 and that otherwise U' is equivalent to U with 
respect to the procedure in Appendix 1 of Reference 3. As a 
consequence of this, the procedure can be extended to include 
station uncertainty by replacing o by [o + ° s ] 2 where ever 
it is used. In this case, o = re where r is the range of an 
initial estimate of an object's position from the assumed station 
position and e is the bearing error (standard deviation) in 
radians of the bearings associated with the station. 
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III. Line of Position Uncertainty and Position Estimates 

The procedure that is described in this section relates 
position estimates that are based on lines of position to line 
of position uncertainty. The procedure is based on a model that 
is defined as follows: Lines of position are straight lines. 

Observed lines of position are parallel to true lines of posi- 
tion. The distance of an observed line of position from a true 
line of position is determined by an independent normally dis- 
tributed random variable with known mean and standard deviation. 
Lines of position are specified in terms of a rectangular coor- 
dinate system with the origin at a reference point as shown in 
Figure 4. For celestial navigation, an appropriate choice for 
the reference point would be the assumed position. Since bear- 
ing lines are lines of position, this model differs from the 
model that is described in Section II only in terminology. 
However, operationally the use of the model that is described 
in this section differs in the way the standard deviation of the 
distance of the line from a true line is determined. The stan- 
dard deviation associated with each line of position must be 
specified. If this can be done, the procedure can be used. 

As an example, suppose that the values are for the first 
line of position and for the second where 0-^ > o ^ and 
that the lines of position intersect at a 90° angle. In this 
example, the minimum area confidence (probability) region is 
an ellipse that is centered on the estimated position. And the 
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position 



o 



1st line of position 



initial estimate 
and final estimate 




2nd line of position 



position 



Figure 4. The geometry when only two lines of position are used 
with the procedure in Appendix 1 of Reference 3. The lines of 
position correspond to bearing lines for observed bearings ^ and 
0£ with respect to the procedure. The procedure determines the 
size and the orientation of elliptical confidence regions of con- 
fidence p, the location of the initial estimate and the location 
of the final estimate. The locations of the estimates correspond 
when only two lines are used. If the lines were obtained from 
sextant observations, the assumed position should be the reference 
position in which case the lines of position would be determined 
by the azimuth anges ot^ and and the distances , ^ and p 2 * 
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estimated position is at the intersection of the bearing lines. 

2 

For a confidence (probability) of containment of 1 - exp(-k /2), 
the major axis of the ellipse is coincident with the first bear- 
ing line and it is of length 2ko^, the minor axis of the ellipse 

is coincident with the second bearing line and it is of length 

2 

2 k 02 and the area of the ellipse is Trk o ^ . 
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IV. A Composite Position Estimate 

The procedure that is described in this section is for 
combining position estimates for an object from independent 
sources. It is based on the following model: The rectanqular 

coordinates of each position estimate are determined by an inde- 
pendent bivariate normal distribution whose covariance matrix 
is known but whose mean vector is not known. The components of 
the mean vector for each of the distributions are the unknown 
coordinates x and y of the object. This model implies that the 
natural logarithm of the likelihood function for a set of n 
estimates can be expressed as follows: 

n i -1 

log L = K - 1/2 Z (x. - x) l. (x, - x) 

i=l 1 1 

where K is a constant, x^ is an estimate vector with components x^ 

and y^, x is the common mean vector with components x and y, 

the unknown coordinates of the object, and Y ■ is the covariance 

2 2 

matrix with elements a.~ , a.- and . The maximum likeli- 

i y lxy 

hood estimates x and y of the unknown coordinates x and y are 
the solutions of the two simultaneous linear equations deter- 
mined by 



9 ( log L) 

3x 



x=x 

/N. 

y=y 



0 and 



9 ( log L) 
3y 



x=x 

y=y 



o 
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The equations can be written as 



A x + B y = D 

B x + C y = E 



and their solution as 



x = 



CD ~ BE 
AC - B 2 



y = 



AE ~ BD 
AC - B 2 



where A = Ea . , B = Zb., C = Ec . , D = E(a.x. + b.y.), 

1 1 1 li i J i 

E = E(b.x. + c.y.), a. = o.~ /d . , b. = - o-~~/d., c. = o-^/d. 

i i i J i ' l ly ' l l lxy i' i ix ' i 



2 2 2 

d. = o . * a.* - and all of the sums are for i from 1 to n 

i ix ly lxy 



Since they are linear combinations of the estimates and 

y^ , the estimates x and y are determined by a bivariate normal 

distribution. Consequently, all that is required to determine 

this distribution is its mean vector with components and 

x y 

2 2 

and its covariance matrix with elements , o~ and The 

x y xy 

mean vector is determined by 



M- = E{(CD - BE) /(AC - B 2 )} = x 

X 

and 

p- =E{(AE - BD ) / (AC - B 2 )} = y 



And, the covariance matrix is determined by 



o- 2 = E{(CD - BE) - E (CD - BE)} 2 / (AC - B 2 ) 2 

X 

= {C 2 (F+I+2H) - 2CB (H+K+M) + B 2 (G+J+2K) }/ (AC - B 2 ) 2 
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a- 2 = E { ( AE - BD) - E (AE - BD)} 2 /(AC - B 2 ) 2 



{A 2 (G+J+2N) - 2AB (H+K+M) + B 2 (F+I+2L) }/ (AC - B 2 ) 2 



0 — = E{[ (CD-BE) -E (CD-BE) ] * [ (AE-BD) -E (AE-BD) ] }/(AC - B 2 ) 2 
xy 

= { (AC+B 2 ) (H+K+M) -CB(F+I+2L)-BA(G+J+2N) }/(AC - B 2 ) 2 



2 2 

where F = Za ■ o . ~ , G 



. u . ^ 
1 IX 



v , 2 2 
Zb • o . - 

1 IX 



H = Za . b . o . - 



1 1 IX 




J 








where all of the sums are for i = 1 to n . 



By using arguments given in Appendix 1 of Reference 3, one 

can show that the axes of the elliptical confidence regions 

associated with x and y are coincident with an x ' y ' -coordinate 

system where the transformation from the xy-coordinate system 

to this system is the coordinate axes rotation through the 

2 2 

angle y defined by tan 2y = 2 (o~ - o~ ) . For a confidence 

^ ? xy' y x 

p, from Reference 3, the minimum area confidence region is an 

2 

ellipse with semi-axes k o~, anc * ^ °G. ' and area wk a~,crg, 

x y ; x y 

where k = [-2 ln(l-p)] 2 , , 

2 2 2 . 2 2 
o~ . = a~ cos Y - 2o-' / ' cos v sin v + sin ) , 

x 1 x ' xy Y 

2 2 2 2 2 
a-', = a~ sin y + 2cs~~ cos y sin y + a A cos y 

y ' x xy y 

and the center of the ellipse is at the point (x,y). In this 

coordinate system, ~ , = 0 . 

x y 
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The above equations can be used to specify a composite posi 

tion estimate in terms of the location, orientation and size of 

an elliptical confidence region which is generally the form in 

which position estimates of the kind that are being considered 

2 2 

here are specified. But, since values of o~ , o'' and 

^ x y xy 

are required for each of the n estimates that is being combined 
a way is needed for determining these values given the orienta- 
tion and size of an elliptical confidence region. A procedure 
to do this when the orientation is given in terms of the direc- 
tion 6 of the major axis and the size is given in terms of the 
lengths SMJ and SMI of the semi-major and semi-minor axes and 
the confidence p is described next. 

By using an xy-coordinate system in which the positive 
y-axis direction is north and the positive x-axis direction is 
east and with the convention 0° < 6 < 180°, the dependence of 
the value of the rotation angle Y and of the order relation 

between 0% and o~ , on the value of the major axis direction 6 
x y 

is indicated by the following table: 

0° < 6 < 45°: y = 6 and o~ , > o~ , 

' y ' x ' 

45° < 6 < 135°: y = 6 - 90° and o~, > o~, 

x ' y ' 

135° < 6 < 180°: y = 6 - 180° and o~ , > o~, 

- ' y ' x 

With an order relation and a value for p, values for o~, and 

X 

o'', can be determined with values for SMJ and SMI. With values 

y 

2 2 

for o~, , o'', and y, values for o~ , o~ and o^ can be found 

x ‘ y x y xy 

from the following equations: 
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o~ 

X 



2 



2 2 2 2 2 
= o~ , cos Y + a~, sin y, 
x ' y ' ' ' 



2 2 2 2 2 
=o^. sin y + o^, cos v 
y x ' 1 y ' 



and 

2 2 

= (o~, - O'" ) sin Y cos y 

xy y ' x ' 

which can be obtained by inverting the equations above for 3 -, 
and , 2 . 

y 

As an example, suppose the data in the following table 
represent three independent position estimates: 





A 

X 


/s 

y 


6 


SMJ 


SMI 


k 


1 st 


- 3.7 


18.1 


59° 


36 


20 


2 


2 nd 


11.8 


8 . 4 


105° 


37 


11 


2 


3rd 


0 


0 


146° 


45 


23 


2 


Here 


, distances 


are in 


nautical 


miles and p = . 


86 in each 


case 


For 


this example with 


values in 


square nautical 


miles : 






°is 2 - 264 - 


58 , 0. - 
ly 


2 = 159.42 


and 0 ~ ~ = 98 . 
lxy 


89 





o 0 - 2 = 321.35, o ~ 2 = 51.15 and o -- - -78 

2 y ' 2 y 2 xy 

0 - 2 = 249.20, o / 2 = 389.30 and o = -173.4 
3x ' 3y 3xy 

These values give the following composite estimate: 
x = -2.46 nautical miles and y = 12.26 nautical miles. 
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6.21 nautical 



For this case, = 16.92 nautical miles, o~ = 

miles and = -55.1 square nautical miles. And, for k = 2: 

xy 

SMJ = 33.85 nautical miles, SMI = 12.42 nautical miles and 
6 = 102°. The composite confidence region and its three com- 
ponent confidence regions are shown in Figure 5. 

As a second example, suppose each position estimate is deter- 
mined by a circular normal distribution. Then o-~ = o.~ = a. 

1 lx ly l 

and = 0 for i = 1 to n. In this case, the composite 

estimate is: 

x = (Ex./o i )/(Z l/o i ) , y = (Zy i /o i )/(Il/o i ) , 

a- 2 = n/(El/o i ) 2 , a- 2 = n/ (Zl/o^) 2 and = 0. 

In this example, since x and y are determined by a circular 
normal distribution, the minimum area confidence regions are 
circles and orientation is not an issue. 
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Figure 5. The ellipses represent the confidence ellipses of 
the first example. The position estimates are at the center 
of the ellipses. The numbers indicate the order of the esti- 
mate in the table on Page 17. 
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